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Starting from the 3D Gross-Pitaevskii equation and using a variational approach, we derive an 
effective ID wave-equation that describes the axial dynamics of a Bose condensate confined in an 
external potential with cylindrical symmetry. The trapping potential is harmonic in the transverse 
direction and generic in the axial one. Our equation, that is a time-dependent non-polynomial 
nonlinear Schrodinger equation (ID NPSE), can be used to model cigar-shaped condensates, whose 
dynamics is essentially ID. We show that ID NPSE gives much more accurate results than all other 
effective equations recently proposed. By using ID NPSE we find analytical solutions for bright and 
dark solitons, which generalize the ones known in the literature. We deduce also an effective 2D non- 
polynomial Schrodinger equation (2D NPSE) that models disc-shaped Bose condensates confined in 
an external trap that is harmonic along the axial direction and generic in the transverse direction. 
In the limiting cases of weak and strong interaction, our approach gives rise to Schrodinger-like 
equations with different polynomial nonlinearities. 
PACS Numbers: 32.80. Pj; 67.65. -f-z; 51.10.-fy 

I. INTRODUCTION 

Bose condensates are nowadays routinely produced by 
many experimental groups all over the world. The ther- 
mal and dynamical properties of Bose condensates have 
been investigated with different atomic species and var- 
ious trap geometries. An interesting theoretical prob- 
lem is the derivation of ID and 2D equations, describ- 
ing cigar-shaped and disc-shaped condensates, respec- 
tively. At zero temperature, a good theoretical tool for 
the study of the dynamics of dilute condensates is the 
time-dependent 3D Gross-Pitaevskii equation . In the 
case of reduced dimensionality, various approaches have 
been adopted to derive effective equations from the 3D 
Gross-Pitaevskii equation (2|-U]. 

In this paper we analyze both cigar-shaped and disc- 
shaped condensates. By using a variational approach, we 
obtain an effective ID time-dependent non-polynomial 
nonlinear Schrodinger equation that describes the axial 
dynamics of a Bose condensate confined in an external 
potential with cylindrical symmetry. We demonstrate 
that our equation exactly reproduces previous findings in 
the limits of weak-coupling and strong-coupling. More- 
over, we show that our variational approach is more ac- 
curate than all other recently proposed procedures in 
the evaluation of both static and dynamical properties 
of the condensate. We also investigate the 2D reduction 
of the 3D Gross-Pitaevskii equation. In this case it is 
not possible to analytically determine an single effective 
2D wave-equation which describes the dynamics of disc- 
shaped condensates. Nevertheless, analytical equations 
can be found in the weakly-interacting limit and in the 
strongly- interacting limit. 



II. EFFECTIVE ID EQUATION 



The 3D Gross-Pitaevskii equation (3D GPE), which 
describes the macroscopic wavefunction ip{r,t) of the 
Bose condensate, is given by 



2m 



U{r)+gN\ij{r,t)\' 



^(r,t) , 
(1) 



where U{r) is the external trapping potential and g — 
as/ra is the scattering amplitude and the s-wave 
scattering length |jl|]. TV is the number of condensed 
bosons and the wave-function is normalized to one. Note 
that the 3D GPE is accurate to describe a condensate of 
dilute bosons only near zero temperature, where thermal 
excitations can be neglected ||] . 

The 3D GPE can be obtained by using the quantum 
least action principle, i.e. 3D GPE is the Euler-Lagrange 
equation of the following action functional 



S = / dtdv ?/'*(r,t) 



dt ^ 2^' 



C/(r) 



V'(r,t) . 



(2) 



We consider an external potential with cylindrical sym- 
metry. In particular we analyze a trapping potential that 
is harmonic in the transverse direction and generic in the 
axial direction: U{r) = ^muj'^{x'^ + y'^) + V (z) . We want 
to minimize the action functional S by choosing an ap- 
propriate trial wavefunction. A natural choice Q| is the 
following 
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ip{r, t) = (t>{x, y, t; a{z, t)) f{z, t) , 



(3) 



where both (p and / are normahzed and is represented 
by a Gaussian: 



e 2<T(i,t)^ 
7ri/V(z,t) ■ 



(4) 



The variational functions a{z, t) and /(z, t) wiU be deter- 
mined by minimizing the action functional after integra- 
tion in the {x, y) plane. The choice of a Gaussian shape 
for the condensate in the transverse direction is well jus- 
tified in the limit of weak interatomic coupling, because 
the exact ground-state of the linear Schrodinger equation 
with harmonic potential is a Gaussian. Moreover, for the 
description of the collective dynamics of Bose-Einstein 
condensates, it has already been shown that the varia- 
tional technique based on Gaussian trial functions leads 
to reliable results even in the large condensate number 
limit [||. 

We assume that the transverse wavefunction (f> is slowly 
varying along the axial direction with respect to the 
transverse direction ||], i.e. ~ V5_</i where = 

+ By inserting the trial wave- function in (2) 

and after spatial integration along x and y variables the 
action functional becomes 
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The Euler-Lagrange equations with respect to /* and a 
read 
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The second Euler-Lagrange equation reduces to an al- 
gebric relation providing a one to one correspondence 
between a and /: ct^ = a']_^l + 2asN\f\'^, where a± — 

\ ^ is the oscillator length in the transverse direction. 
One sees that a depends implicitly on z and t because 
of the space and time dependence of Inserting this 
result in the first equation one finally obtains 
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This equation is the main result of our paper. It is a time- 
dependent non-polynomial nonlinear Schrodinger equa- 
tion (ID NPSE). 

We observe that from ID NPSE in certain limit- 
ing cases one recovers familiar results. In the weakly - 
interacting limit as-/V|/p << 1 one has 
previous equation reduces to 



a]_ and the 



2m dz'^ 



V - 



gN 
27raT 



f 



(9) 



where the additive constant hLUj_ has been omitted be- 
cause it does not affect the dynamics. This equation 
is a ID Gross-Pitaevskii equation. The nonlinear co- 
efficient g' of this ID GPE can be thus obtained from 
the nonlinear coefhcient g of the 3D GPE by setting 
g' = g/{2na\). This ansatz has been already used by 
various authors, for example in Ref. Note that the 
limit CsA^I/P << 1 is precisely the regime where the 
healing length is larger than a. In this regime the cigar- 
shaped condensate is quasi- ID, as shown in a recent ex- 
periment HJ. It is well known that in ID and at ultra- 
low densities (ajA^j/p << a1/a\ « 1) an interacting 
Bose gas becomes a Tonks gas, i.e. a gas of spinless 
Fermions Such a transition cannot be described by 
the cubic 3D GPE equation and therefore by ID NPSE 
because in the Tonks regime the inter-atomic interaction 
cannot simply be approximated by a zero-range pseudo- 
potential in mean field approximation pO[ |. Instead, in 
the strongly- interacting high density limit asiV|/p >> 1 
(but iV|-0pas << 1 to satisfy the diluteness condition) 
one finds = y/2a\aj^ N^l^\f\ and the ID NPSE be- 
comes 
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In this limit, and in the stationary case, the kinetic term 
can be neglected (Thomas- Fermi approximation) and one 
finds the following analytical formula for the axial density 
profile 
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where /i' is the chemical potential, fixed by the normal- 
ization condition. It is important to stress that this ID 
Thomas-Fcrmi density profile is quadratic in the term 
fi' — V{z). The same quadratic dependence is obtained 
starting from the Thomas-Fermi approximation of the 3D 
stationary GPE, i.e. neglecting the spatial derivatives in 
Eq. (1), and then integrating along x and y variables. In 
this way one finds a formula that differs from (11) only 
for the numerical factor which is 1/4 instead of 2/9. 
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To test the accuracy of the full ID NPSE, Eq. (8), 
and compare it with other procedures proposed in the 
last few years, we numerically investigate the simple 
case of harmonic trapping also in the axial direction: 
V{z) = imw^z^. In this case, the ID NPSE can be 

written in scaled units: z in units of a, = i/— the 
oscillator length in the axial direction, and t in units 
oi 1/ujz. In order to assess the accuracy of the various 
ID approximations we have also solved numerically the 
3D Gross-Pitaevskii equation (3D GPE), given by Eq. 
(1) with imaginary time [0. Note that the numerical 
solution of the ID NPSE is not more time consuming 
than the solution of the standard ID GPE. In Fig. 1 we 
plot the normalized density profile p{z) = of the 

ground-state wave-function of a cigar-shaped condensate 
confined in a trap with an aspect ratio az/a± — 10. 
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FIG. 1. Normalized density profile p{z) — along the 

axial direction z for the cigar shaped trap. Number of Bosons: 
N = 10'* and trap anisotropy: u)±/u)z = 10. Four differ- 
ent procedures: 3D GPE (solid line), ID GPE (dashed line), 
CGPE (long-dashed line) and ID NPSE (dotted line). From 
top to bottom: as/flz = 10~*, as/a^ — 10~^, as/a^ — 10^"^. 
Length z in units of and density in units of . 

We compare the results of four different procedures. 
The first procedure is the "exact" one, i.e. the solution 
of the 3D GPE. The second procedure is the numerical 
solution of the ID Gross-Pitaevskii equation (ID GPE) 
given by Eq. (9) with an imaginary time. The third 
procedure is that proposed by Chiofalo and Tosi Q. In 
this case the nonlinear term of the ID Gross-Pitaevskii 
equation is found by imposing that the ID wave-function 
has the same chemical potential of the 3D one (CGPE). 



The fourth and last procedure is the numerical solution of 
our non-polynomial nonlinear Schrodinger equation (ID 
NPSE), i.e. Eq. (8) with imaginary time. As shown in 
Fig. 1, the ID NPSE results are always very close to 
the "exact" ones and much better of the other approx- 
imations. Moreover, the CGPE procedure gives better 
results than the ID GPE for large values of the scatter- 
ing length where Eq. (9) is not reliable but in any case 
the ID NPSE is superior. 

The very good performance of Eq. (8) is not limited 
to the ground-state. We have investigated the dynam- 
ics of the condensate by taking the previously calculated 
ground-state wavefunctions but changing the harmonic 
trap in the axial direction: from V{z) = ^mujlz'^ to 
V{z) = ^mcOzz'^. In this way the condensate shows large 
collective shape oscillations along the z axis. 
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FIG. 2. Squared amplitude < 2^ > as a function of time 
t. Number of Bosons: A*' = IC* and trap anisotropy: 
ij±/ujz — 10. Four different procedures: 3D GPE (solid 
line), ID GPE (dashed line), CGPE (long-dashed line) and 
ID NPSE (dotted line). Scattering length: a^/a^ = 10"^ 
Length z in units of a^, density in units of and time t in 
units of 1/uJz. 

In Fig. 2 we plot the time evolution of the squared 
amplitude < > of the condensate in the axial direc- 
tion, given by < >= J dzz"^ p{z,t), where p{z,t) = 
J dxdy\il;{r,t)\'^ in the case of the 3D GPE and p{z,t) — 
|/(z,i)P in the other cases. Apart the better evaluation 
of the amplitude that is a consequence of the better eval- 
uation of the ground-state wave- function, one sees that 
our ID NPSE reproduces quite well also the "exact" sinu- 
soidal behavior of the collective oscillation. For example, 
after four oscillations the relative error in the determina- 
tion of the instant t of minimum radius < > is about 
1%. The other two procedures (CGPE and ID GPE) 
clearly show a frequency delay with respect to the 3D 
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GPE and ID NPSE results: after four oscillations their 
relative errors, with respect to the 3D GPE result, are 
about 5% and 6% respectively. 

We remind that the ID NPSE has been obtained by 
using a factorization of the 3D wavefunction with a vari- 
ational ansatz for the transverse part of the wavefunction 
and neglecting the term ^0*§^- Under the condition 
of the present computations the last assumption is fully 
justified because we have numerically verified that the 
ratio between the neglected term and the total energy 
ranges form 10"'^ in the weak-coupling limit to 10~^ in 
the strong-coupling limit. 
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FIG. 3. Four frames of the axial density profile p{z) of the 
Bose condensate tunneling through the Gaussian barrier (Eq. 
(3)). Start-up kinetic energy per particle of the condensate: 
Eo = 10. Gaussian barrier parameters: Vb = 10 and E = 1. 
Interaction strength: Na^/a^ = 10. Comparison among 3D 
GPE (solid line), ID NPSE (dotted line) and ID GPE (dashed 
line). Length in units of — {h/muJ:,)^^'^ , density in units of 
, time in units of and energy in units of hujz. 

To test the accuracy of ID NPSE in the description 
of the dynamics of a cigar-shaped Bose condensate in 
more complex problem, we investigate the scattering and 
tunneling of the condensate on a Gaussian barrier. The 
initial wave function of the condensate is found by solv- 
ing the equations with imaginary time and imposing a 
harmonic trapping potential also in the horizontal axial 
direction: 



V{z) = Knujl{z~ Zq) 



(12) 



To have a cigar-shaped condensate we choose lo^Ilo^ = 
10. We set zq = 20, where zq is written in units of the 
harmonic length = {h/mWzY^^ . For t > the trap in 
the axial direction is switched off and a Gaussian energy 
barrier is inserted at z = 0. The potential barrier is given 
by 



Viz) = Vo e-^''/^' 



(13) 



where Vq is the height of the potential barrier and S its 
width. The condensate is moved towards the barrier by 
adding an initial momentum po in the axial direction: 



f{z,Q)~. f{z,{))e'^P^^'^ 



(14) 



In the case of a non-interacting condensate, i.e. as = 0, 
the total energy per particle is practically the start-up 
kinetic energy i?o = pQ/(2m), and the transverse energy 
TiLj±^ does not affect the dynamics. 

As shown in Fig. 3, where we plot the density profile 
of the condensate along the symmetry axis at different 
instants, a fraction of the condensate tunnels the barrier 
while the rest is reflected. As expected, one sees also the 
interference between the incident wave function and the 
reflected wave function. Fig. 3 shows that the "exact" 
axial density profile obtained with the 3D GPE and that 
of the ID NPSE are always quite close, also during the 
impact and tunneling time. These results suggest that 
ID NPSE is very adequate also in the description of tun- 
neling phenomena. 

The cigar-shaped configuration of the condensate is 
useful to study topological objects, like bright and dark 
solitons. By using ID NPSE one finds out solitonic so- 
lutions which generalize what has been previously found 
with ID GPE |,|l|. Dark solitons (a^ > 0) of Bose con- 
densed atoms have been experimentally observed p3| , 
while bright solitons (os < 0) are more elusive due to the 
collapse of the condensate with a large number of atoms. 

Let us first consider bright solitons. Starting from our 
ID NPSE, setting V{z) = 0, scaling z in units of a± and 



t in units of < 



with the position 



f{z,t)^<^{z-vt)e' 
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we find 




/i* . (16) 



where Q = z — vt and 7 = \as\N / a±. This is a Newto- 
nian second-order differential equation and its constant 
of motion is given by 
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The boundary condition $ — > for C — > cx) implies that 
E — Q. Then, by quadratures, one obtains the bright- 
soliton solution written in implicit form 
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Moreover, by imposing the normalization condition one 
also finds 



(1 



7 = 



(19) 



The normalization relates the chemical potential to 
the coupling constant 7, while the velocity v of the 
bright soliton remains arbitrary. In the weak-coupling 
limit (7$'^ << 1), the normalization condition gives 
1 — 7^/2 and the bright-sohton solution reads 



*(C) = A/^sec/i[7C] 



(20) 



The above solution is the text-book bright soliton of the 
ID non-linear cubic Schrodinger equation (ID GPE). As 
shown in Figure 4, for ^/2/3 < 7 < 2/3 there are two 
values for the chemical potential but we have numeri- 
cally verified that lower one corresponds to an unstable 
solitonic solution. For 7 > 2/3 there are no solitary- wave 
solutions due to the collapse of the condensate. 



the weak-coupling and in the strong-coupling limit. In 
the weak-coupling limits one finds 



(21) 



where $0 = 7^/^/V2 and 7 = \as\l\N/a^, with AiV 
the number of missing Bosons due to the hole in the 
condensate produced by the dark soliton. Eq. (20) is 
the familiar formula of stationary dark soliton for the ID 
nonlinear Schrodinger equation (ID GPE), while in the 
strong-coupling limit one has 



(22) 



where a = V37^/'*V^/2^/^, b = arcth[l / \f?,] and $0 = 
71/6/(72(73-2/3)2/3). 

The study of the stability of our solitary-wave solu- 
tions is left to a future work. 3D numerical calculations 
[Q suggest that Bose condensed dark solitons are stable 
for sufhciently small numbers of atoms or large trans- 
verse confinement. Moreover, a recent experiment p^ ] 
has shown that dark solitons, created in a spherical Bose 
condensate, decay into vortex rings. 
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FIG. 4. On the top: chemical potential ^ of the bright soli- 
ton as a function of the interaction strength 7 = \as\N/a±. 
ID NPSE (full line) and ID GPE (dashed line). On the bot- 
tom: axial density profile p((") = $^(C) of the bright soliton 
of Bose condensed atoms for three values of the interaction 
strength 7 obtained with ID NPSE. 

In the case of static dark solitons, the boundary con- 
ditions are: <i>(0) = and <i> — > $0 for 2; — > 00. Starting 
from the ID NPSE the analytical formula of the dark 
soliton is quite intricate because it involves elliptic inte- 
grals. Nevertheless, simple expressions can be found in 



III. EFFECTIVE 2D EQUATION 

In this section we derive the 2D reduction of the 3D 
GPE equation by using again a variational approach. In 
this case we take a trapping potential that is harmonic in 
the axial direction and generic in the transverse direction: 
^(r) = W{x^y) -\- ^mLolz^ . As trial wave-function we 
take the following 



V^(r, t) = (f){x, y, t) f{z, t; r]{x, y, t)) 



with 



f{z,t;r]{x,y,t)) = 



7ri/4j^(a;^y,f)i/2 
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where r]{x, y, t) is a variational function that describes the 
width of the condensate in the axial direction. We assume 
that the axial wavefunction / is slowly varying along the 
transverse direction with respect to the axial direction, 
i.e. V^/ ~ By inserting the trial wave-function 

in (2) and after spatial integration along z variable, one 
finds the following Euler-Lagrange equations: 
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2(27r)^/^ 

Note that the second Euler-Lagrange equation can be 
written as rf ~ 2 (27r)^^^ of a^A^l^pTy — of = 0, and, con- 
trary to the case of ID reduction, it does not have a el- 
egant analytical solution but it can be easily solved . 
We name the 2D non-polynomial Schrodinger equation 
(25), with the condition (26), 2D NPSE. In Fig. 5 we 
compare the radial density profile p{r) of the Bose con- 
densate obtained by solving 2D NPSE with the exact one 
obtained by solving the 3D GPE. For simplicity, we use 
a harmonic trapping potential W{x, y) — ^muj'j^{x'^ + y'^) 
also in the transverse radial direction. Fig. 5 shows that 
2D NPSE is quite reliable in describing the ground-state 
of disc-shaped Bose condensate. For example, the rela- 
tive error in the calculation of the density at the orgin 
ranges from 0.7% to 5% by increasing asja^. 
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FIG. 5. Radial density profile p(r) of the Bose condensate 
for a disc shaped trap. Number of Bosons: A'^ = 10** and 
trap anisotropy: uj±/ujz ~ 1/10. 3D GPE (solid line) and 
2D NPSE (dotted line). From top to bottom: Us/a^ = 10"", 
ds/ciz = 10~^, as/az = 10~^. Length r — + y'^ in units 
of and density in units aj^. 

We can analyze the limits of weak and strong inter- 
action of 2D NPSE. In the weakly-interacting case one 
finds J] ~ az and the equation of motion becomes 
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where the constant hLo^ has been omitted because it does 
not affect the dynamics. This equation is a 2D Gross- 
Pitaevskii equation. The nonlinear coefficient g" of this 
2D GPE can be thus obtained from the nonlinear coef- 
ficient g of the 3D GPE by setting g" = g / {{2TTf/^az). 



Eq. (27) describes a disc-shaped Bose condensate in a 
quasi-2D configuration. In fact, the weakly-interacting 
limit corresponds to a condensate with a chemical po- 
tential lower than fiuj^ [^j. In the strongly-interacting 
case one has instead 77 = V^ii^at^^ aV^ N^/^\(j)\'^^'^ , and 
the resulting nonlinear Schodinger equation is 
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Note that in this limit, and in the stationary case, the 
kinetic term can be neglected (Thomas-Fermi approxi- 
mation) and one finds the following analytical formula 
for the transverse density profile 
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{^Ji" -W{x,y)) 
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(29) 



where ^" is the chemical potential, fixed by the nor- 
malization condition. The same dependence is obtained 
starting from the Thomas-Fermi approximation of the 
stationary 3D GPE, i.e. neglecting the spatial deriva- 
tives in Eq. (1), and then integrating along z variable. 
In this way one finds a formula that differs from the pre- 
vious one only for the numerical factor which is \/2/(37r) 
instead of l/(3\/37r). 



CONCLUSIONS 

We have found that the ID non-polynomial nonlinear 
Schrodinger equation we have obtained by using a Gaus- 
sian variational ansatz from the 3D Gross-Pitaevskii ac- 
tion functional is quite reliable in describing the ground 
state and axial collective oscillations of cigar-shaped con- 
densates. We have also tested the accuracy of the non- 
polynomial nonlinear equation in the case of a Bose 
condensate scattering and tunneling on a Gaussian bar- 
rier. The agreement with the results of the 3D Gross- 
Pitaevskii equation is very good for both ground-state 
and dynamics of the condensate. This equation will be 
useful for detailed numerical analysis of the dynamics 
of cigar-shaped condensates, particularly when the local 
density may undergo sudden and large variations. In 
fact, the accurate mapping we have provided, allows to 
maintain a very good spatial resolution with modest com- 
putational effort even if the breakdown of the weak in- 
teraction condition during time evolution prevents the 
use of standard ID GP equation. This case is often en- 
countered in the study of reflection and tunneling events 
and in the propagation of solitary waves, which are cur- 
rently being experimentally investigated. By using the 
ID non-polynomial Schrodinger equation we have ob- 
tained analytical formulas for bright and dark solitons 
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which generahze the ones known in the hterature for 
the ID non-Hnear cubic Schrodinger equation. Finally, 
we have deduced effective 2D equations describing the 
axial dynamics of disc-shaped condensates. From the 
3D Gross-Pitaevskii action functional and using another 
Gaussian variational ansatz we have derived the effec- 
tive 2D nonlinear Schrodinger equation, which again has 
a non-polynomial structure. This structure simplifies in 
the weak and strong interacting regime. 
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